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We present the constraints on isocurvature fluctuations for a flat universe implied by the Boomerang 
and Maxima-1 data on the anisotropy of the cosmic microwave background. Because the new data 
defines the shape of the angular power spectrum in the region of the first acoustic peaks much more 
clearly than earlier data, even a tilted pure isocurvature model is now ruled out. However, a mixed 
model with a sizable isocurvature contribution remains allowed. We consider primordial fluctuations 
with different spectral indices for the adiabatic and isocurvature perturbations, and find that the 
95% C.L. upper limit to the isocurvature contribution to the low multipoles is a < 0.63. The upper 
limit to the contribution in the I ~ 200 region is Q200 < 0.13. 
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I. INTRODUCTION 

Two recent balloon-borne experiments, Boomerang Jl]] 
and Maxima-1 ||, have lent considerable support to 
the inflationary paradigm. The Boomerang experiment, 
which was flown over Antarctica for 10 days in 1999 and 
which measured the temperature fluctuations over 1% 
of the microwave sky, was the first to provide firm evi- 
dence for the first (and possibly second) acoustic peak. 
The angular power spectrum, obtained from a prelimi- 
nary analysis of the Boomerang data, is compatible with 
approximate scale invariance predicted by cosmic infla- 
tion. The first acoustic peak is at the multipole I ~ 200, 
implying a flat universe with f2 = f2 m + ~ 1 in the 
case of adiabatic fluctuations Q. The Boomerang re- 
sults have been confirmed by Maxima-1 which ob- 
served a 124 square degree patch of the sky. Boomerang 
and Maxima-1 are the first precision measurements of 
the cosmic microwave background (CMB) temperature 
fluctuations at multipoles Z > 30. They will be followed 
by the second generation satellite experiments Microwave 
Anisotropy Probe (MAP) j§ and Planck ||. 

The simplest one-field inflation model typically pre- 
dicts a near scale invariant power spectrum with Gaus- 
sian fluctuations. The perturbations generated by quan- 
tum fluctuations are adiabatic, with the number density 
proportional to entropy density so that 8{n/s) = 0. This 
is so because the quantum fluctuations of the inflaton 
field are directly reflected in perturbations of the infla- 
ton energy density. 

More generally, in addition to adiabatic fluctuations, 
there is also another fundamental fluctuation mode. It 
does not perturb the total energy density on comoving 
hypersurfaces (orthogonal to fluid flow), so that 5p = 0. 



Since there is then no curvature perturbation of the co- 
moving hypersurfaces, it is called an isocurvature fluctu- 
ation m?]] . Such fluctuations can arise in particle physics 
models where, during inflation, there appear also other 
low mass fields in addition to the inflaton. For instance, 
for a complex scalar field with a U(l) symmetry the fluc- 
tuations of the phase do not affect the energy density 
but would show up as perturbations in the conserved 
charge associated with the symmetry. This is effectively 
the case in the minimally supersymmetric model, and 
its extensions, which have several flat directions in the 
potential by virtue of gauge invariance and supersymme- 
try. The complex fields along these flat directions, called 
Affleck-Dine fields 0, will be subject to quantum fluctu- 
ations similar to the inflaton. As the Affleck-Dine fields 
typically carry global charges such as baryon number, 
in addition to adiabatic perturbations, there will also be 
isocurvature perturbations || . The Affleck-Dine conden- 
sate itself is not stable but fragments into non-topological 
solitons 1 10 1, or Q-balls. These may be absolutely stable, 
in which case there will be isocurvature fluctuations in 
the baryon number, or unstable but long-lived, in which 
case isocurvature fluctuation can be imprinted also on 
the cold dark matter (CDM) particles pip . Other par- 
ticle physics motivated sources for isocurvature fluctua- 
tions are axion models fl^| or models of inflation with 
more than one field fll3| . 

For isocurvature perturbations the overdensities in 
a given particle species are balanced by perturbations 
in other particle species, such as radiation. At the 
last scattering surface (LSS) the compensation for the 
isocurvature perturbations can be maintained only for 
scales larger than the horizon, effectively generating ex- 
tra power to photon perturbations at small multipoles 
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I. As a consequence, the CMB angular power spectrum, 
Ci , of isocurvature perturbations differs a great deal from 
adiabatic perturbations, and a purely isocurvature CDM 
perturbation with a scale free spectrum is clearly ruled 
out H] on the basis of the Cosmic Background Explorer 
(COBE) normalization jl5],[l6| and erg, the observed am- 
plitude of the rms mass fluctuations in an 8/i _1 Mpc 
sphere. 

To match as with the COBE normalization requires a 
large "blue" tilt n; so w 2.2. Such a tilted isocurvature 
model can be obtained from a reasonable inflation model 
JItJ and agrees with most observational data. (The main 
difficulty before Boomerang was that it has more tilt at 
low multipoles than what COBE saw.) Thus it could be 
argued that a pure isocurvature model was still accept- 
able @. 

A mixed perturbation, with a small isocurvature com- 
ponent is certainly allowed and indeed motivated by par- 
ticle physics, although the forthcoming satellite experi- 
ments are expected to be able to constrain the isocurva- 
ture amplitude with a high precision [ pj|j20| . 

The purpose of the present paper is to find out the 
constraints on the isocurvature component as implied by 
the Boomerang and Maxima- 1 data. The isocurvature 
perturbation is highly degenerate with the tensor per- 
turbation, and for their separation one probably has to 
wait for the polarization data from Planck |l9) . Here we 
will make the simplifying assumption that the tensor con- 
tribution is negligible. This will be the case in so-called 
small- field inflation models 21 2^]. For large-field infla- 
tion, the tensor/scalar ratio is r = Cj/Cf d < 7(1 — n a a), 
so that the tensor contribution disappears for a scale-free 
spectrum, n a d ~ 1. As we shall see, such an assumption 
is consistent with the Boomerang and Maxima- 1 data. 

An isocurvature contribution is also somewhat degen- 
erate with reionization. As we are looking for the up- 
per limit to an isocurvature contribution, we make the 
simplifying assumption that reionization has a negligible 
effect. 



II. ADIABATIC VS ISOCURVATURE 
PERTURBATIONS 

Adiabatic perturbations are characterized by 

the quantity 1Z, which is related to the curvature pertur- 
bation of the comoving hypersurface, 



$c — (3/4)(5 7 , where n c is the number density of CDM, 
s 7 is the entropy density associated with photons, and 
8 C and 5 1 are the relative overdensities in the CDM and 
photon energy densities. During radiation domination 
S ~ 5 C = 6n c /n c . Outside the horizon, S does not change 
with time. 

The terms "adiabatic" and "isocurvature" refer to "ini- 
tial conditions" specified at a time when the universe is 
radiation-dominated and all relevant scales are well be- 
yond the horizon. During this era S and TZ are time- 
independent. An adiabatic perturbation is then defined 
as one for which S = 0, and an isocurvature perturba- 
tion as one for which TZ = 0. A general perturbation is a 
superposition of an adiabatic and an isocurvature pertur- 
bation (in fact, there can be several kinds of isocurvature 
perturbations |Q as there are several kinds of matter; we 
consider here a CDM isocurvature perturbation). 

If the perturbations are Gaussian, their statistical 
properties are fully described by their power spectra 
P n (k) = {\TZ k \ 2 ) and P s (fc) = (\S k \ 2 ). The spectral in- 
dex is defined by n(k) = d\nP(k)/dk + 4. If the indices 
are scale independent, we can write 



P n (k) = Ak n ^~ 4 
P s {k) = Bk n "°~ 4 



(2) 



where n ac j and nj SO are the spectral indices with n ac i = 1 
and ni so = 1 corresponding to a completely scale free 
spectrum. (For n a d, this definition is standard, for 7ij so , 
other definitions in use set the index for a scale- free spec- 
trum to rii so — or ni so = —3.) 

While the perturbations are outside the horizon, S re- 
mains time independent, but an isocurvature perturba- 
tion develops a nonzero TZ = when the universe be- 
comes matter dominated p2| . Thus, during last scatter- 
ing a curvature perturbation is already present. During 
matter domination, the density perturbation is related to 
the curvature perturbation by 



5\aH 



TZi 



(3) 



while outside the horizon (k <C aH). 

For low multipoles, / 200, the corresponding scales 
are well outside the horizon during last scattering and the 
CMB anisotropy has a simple relation to the initial per- 
turbations. For adiabatic perturbations, the anisotropy 
comes from the Sachs- Wolfe-effect, 
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where a is the scale factor, k is the comoving wave num- 
ber, and i?( 3 ) is the scalar curvature of the hypersurface 
(for 0=1, the unperturbed hypersurface is flat). For 
adiabatic perturbations 7Z is independent of time while 
outside the horizon. 

Isocurvature perturbations P ]25| , |2"4j | are characterized 
by the entropy perturbation S = (S(n c /s 7 )/(rt c /s 7 ) = 



5T 1 



(4) 



whereas for isocurvature perturbations there is also a di- 
rect contribution from the entropy perturbation, which 
is 5 times larger, 

T 5 3 15 5 y ' 

(This factor 6 in ST/T becomes a factor 36 in the C/.) 
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For higher multipoles the perturbations have evolved 
by the time of last scattering and the perturbations con- 
tribute to the anisotropy by other mechanisms also. The 
full angular power spectra from given initial perturba- 
tions need to be calculated by computer codes such as 
CMBFAST & 



III. COSMOLOGICAL MODELS 

We allow for different spectral indices ( "tilts" ) n a d and 
rij so for adiabatic and isocurvature perturbations, and 
define 



CMB data becomes available, one may need to seriously 
consider non-standard nucleosynthesis scenarios |3C[ |. 

We first used a coarse grid in the range (following [p6|) 
n A = 0-0.8, uj h = 0.003-0.130, uj c = 0.02-0.80. We then 
refined it to focus on the relevant region; for adiabatic 
and mixed models: Q A = 0.48-0.78, uj b = 0.018-0.042, 
lu c = 0.13-0.21, with about a dozen values for each pa- 
rameter; and for pure isocurvature models: = 0.25- 
0.75, uj b = 0.004-0.028, uj c = 0.10-0.40. For the spectral 
indices we used a denser grid with step 0.02; in the range 
"ad = 0.5-2.0, rij so = 0.5-2.5 for adiabatic and mixed 
models; and rii so = 1.5-2.7 for pure isocurvature models. 



a 



or 



Cf 
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to indicate the relative contribution of the isocurvature 
fluctuations. Here the C 2 so and C| d are the expectation 
values of the isocurvature and adiabatic contributions to 
the quadrupole (C 2 = Cf° + Cf). 

Note that different authors @IM|,|§ US6 different 
definitions for this parameter. Since our definition fixes 
the isocurvature contribution at the low end of the mul- 
tipole range, the relative isocurvature contribution at 
high multipoles can be quite different depending on rii so 
and n a d- Since isocurvature perturbations produce much 
more power at low multipoles, the isocurvature contri- 
bution to high multipoles is much smaller than a for 

"iso ^ "ad ■ 

Motivated by inflation, we restrict ourselves to mod- 
els with Q = 1. We also take r — 0, assume r = 
for the optical depth due to reionization, and consider 
a seven-parameter space of cosmological models. These 
parameters are Qa = 1 — £l m , = £l b h 2 , uj c = O c /i 2 , the 
spectral indices n a d and "i so , and the amplitudes A and 
B of the power spectra [Eq. (0)]. The last two translate 
into the isocurvature contribution parameter a and an 
overall normalization. The Hubble constant is given by 
h = a/(w c + w&)/(l - Cl\). 

In the present paper we assume that the adiabatic 
and isocurvature perturbations are uncorrelated and that 
their spectral indices are constant over the cosmological 
range of scales. This is not necessarily the situation in 
all particle physics models, but the nature of correlation 
between isocurvature and adiabatic amplitudes is very 
much model dependent. In general, such correlations 
would result in bounds more stringent than presented 
here. 

We use the COBE p5| , Boomerang [jjj, and Maxima-1 
H data and find the best-fit models which satisfy h = 
0.45-0.90 and a 8 n° m 56 = 0.43-0.70 (top-hat priors). We 
do not use any prior for the baryon density u>b. 

As noted in p7,p8fl, especially the Boomerang data fa- 
vors LOb ~ 0.03, which is above the range allowed by stan- 
dard big-bang nucleosynthesis, uib — 0.006-0.023 p9[ ]. If 
this situation persists in the future, when more accurate 



IV. RESULTS 

The best-fit adiabatic model (model 1, see Table I) is 
close to scale free, n a d = 0.96, and the fit has \ 2 = 23.0 
for 30 data points and 5 parameters. (We did not attempt 
to refine the parameter values between our grid p oints, 
as this is not the focus of this paper. See p7| , p8| j3]| for 
cosmological parameter estimation in adiabatic models 
for the Boomerang and Maxima-1 data.) 

The best-fit isocurvature model (model 2) has rii so = 
2.10 and x 2 = 126.4, for 5 parameters. Clearly a pure 
isocurvature model is ruled out. 

The best fit mixed (adiabatic+isocurvature) model 
(model 3) has n a d = 0.98, n iso = 2.26, and a = 0.008. 
The fit has \ 2 — 22.5, for 7 parameters. The small im- 
provement in the fit from allowing an isocurvature con- 
tribution indicates that the data does not suggest the 
presence of such a contribution (but neither rules it out). 
The big blue tilt ni so = 2.26 in the best-fit isocurvature 
contribution is not significant, since the fit was almost as 
good for any n\ ao between 0.5 and 2.5. The isocurvature 
contribution is so small that it does not affect the angu- 
lar power spectrum very much, and therefore also its tilt 
is not constrained. The small shifts from changing n; so 
can be approximated by changing the other parameters, 
within the accuracy of the present data. 

For large ni SO the smallness of the parameter a is some- 
what misleading, as we defined it as the contribution to 
the lowest multipole. To represent the isocurvature con- 
tribution in the region of the first peak, let us define 
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^-<ad 
°200 



(7) 



where for the purpose of this paper C200 is the average 
1(1 + l)Cj in the 4th Boomerang bin I = 176-225. For 
the above best-fit model, 0200 = 0.04. Similarly, the 
isocurvature contribution to cr| is 8% in this model. The 
models with smaller n iso that are almost as good fits nat- 
urally have larger a (see Fig. ||) but smaller 0200 • Ac- 
cordingly, the best-fit scale- free ("i so ,"ad) = (1,1) mixed 
model (model 4) is almost as good a fit, % 2 = 23.0 for 
only 5 parameters, but it has a larger a = 0.09 and a 
smaller a 2 oo = 0.002. 
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Model 


x 2 








h 


"-ad 




a 


Q200 


(1) Best-fit adiabatic 


23.0 


0.60 


0.028 


0.19 


0.74 


0.96 








(2) Best-fit isocurvature 


126.4 


0.65 


0.012 


0.26 


0.88 




2.10 






(3) Best-fit mixed 


22.5 


0.68 


0.030 


0.20 


0.85 


0.98 


2.26 


0.008 


0.04 


(4) Best-fit scale-free 


23.0 


0.66 


0.030 


0.19 


0.80 


1 


1 


0.09 


0.002 


(5) Large a (Fig. g) 


26.4 


0.72 


0.036 


0.19 


0.90 


1.20 


0.90 


0.63 


0.012 


(6) Large a?200 (rig- 
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(7) Shifted, best-fit adiabatic 


21.6 


0.68 


0.028 


0.15 


0.75 


0.96 








(8) Shifted, best-fit isocurvature 


116.3 


0.70 


0.010 


0.22 


0.88 




2.18 






(9) Shifted, best-fit mixed 


20.9 


0.72 


0.030 


0.15 


0.80 


1.00 


1.78 


0.05 


0.03 


(10) Shifted, best-fit scale-free 


21.7 


0.72 


0.030 


0.15 


0.80 


1 


1 


0.08 


0.0016 


TABLE I. Cosmological models discussed 


in the text and their parameter 


values. 


"Shifted" 


means 


that the fit 


is done after 



shifting the Boomerang data up and the Maxima- 1 data down by their calibration uncertainty. 



For the 2-er (A^ 2 = 4) upper limit to a we find a < 
0.63 and «200 < 0.13. Within scale-free models the upper 
limits are a < 0.34 and 0200 < 0.010. 

In mixed models with a large rii so , the isocurvature 
contribution to <j% can be significant. The 2-er upper 
limit to the isocurvature contribution to a\ is 35%, but 
within scale-free models only 0.3%. 

In Fig. [I] we show the allowed region in the (ni so ,n a d)- 
plane, and the best-fit values of a for each (ni so ,n a( j). 
In Fig. H we show the allowed region (A% 2 < 3.5) in 
the (ni so ,n ac j,Q:)-space as contours of maximum a on the 
(ni so ,n a( j) plane. This indicates the magnitude of the 
isocurvature contribution which remains allowed by the 
present data for different spectral indices. 
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FIG. 1. The 68% (white), 95% (light gray), and 99.7% 
(medium gray) confidence level regions (Ax 2 = 2.3, 6.2, and 
11.8) on the (riiso,n a d)-plane, and the best-fit values of a for 
each (ni so ,n a d). The best-fit model (model 3) is marked with 
an asterisk (*). 

We see that a blue tilt, n a d > 1, in the adiabatic 
spectrum favors a larger isocurvature contribution, since 
the low-multipole power from the isocurvature perturba- 
tions compensates for this tilt. In fact, the presence of 
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FIG. 2. The 68% confidence level region (A X 2 < 3.5) in 
the (ni so ,n a d,a)-space represented by contours of maximum a. 
We do not show minimum a, except that we show the region 
(between the dotted lines) where the minimum is a — 0. 

an isocurvature contribution makes a larger ro a d accept- 
able than in the pure adiabatic case. A large blue tilt, 
nj SO 3> 1, in the isocurvature spectrum makes it less help- 
ful for this purpose, and the upper limit to a goes down 
again. 

Since a measures the contribution to the lowest mul- 
tipoles, the contribution to high multipoles for a given 
a is larger for large n- lso . Thus, although the limits to 
a are tighter for large ni so , this does not mean that the 
contribution to high multipoles would have to be less. In 
fact, the largest allowed CV200; above 0.1, are for spectral 
indices ni so > 1.8 (see Fig. ||). 

In Fig. ^ we show the angular power spectra for some 
of the models discussed above, together with the COBE, 
Boomerang, and Maxima- 1 data. Figures g and ^ show 
separately the adiabatic and isocurvature contributions 
in models 5 and 6, which have large isocurvature con- 
tributions (large a and large c^ooi respectively) but are 
allowed by the present CMB data. 



4 



1.3- 




0.8 i i i i i i i i i i 

0.6 0.8 1 1.2 1.4 1.6 1.8 2 2.2 2.4 

n. 

ISO 

FIG. 3. Same as Fig. g, but for o?200- 

We see how the isocurvature contribution slightly mod- 
ifies the shape of the angular power spectrum. In model 
6 it fills the minima between the peaks in the adiabatic 
contribution. An isocurvature contribution with a large 
blue tilt could be the reason for the lack of prominence 
of the 2nd peak in the Boomerang data. 

The Boomerang and Maxima teams estimate a cal- 
ibration uncertainty of 10% and 4%, respectively, in 
their measurements. This uncertainty is not included 
in the above results. To study its effect, we multiplied 
the Boomerang measurements by 1.1 and divided the 
Maxima- 1 measurements by 1.04 2 . This brings the mea- 
surements closer to each other, since the Maxima-1 points 
tend to lie above the Boomerang points, with the excep- 
tion of the Z e ff = 147 point. The best-fit mixed model 
becomes n ad — 1.00, 7i iso = 1.78, and a = 0.05 (model 
9), with an improved fit, \ 2 — 20.9. The upper limit 
to the isocurvature contribution is tightened to a < 0.56 
(«200 < 0.13 stays the same), but the results do not 
change qualitatively. 

The data is consistent with a purely adiabatic pertur- 
bation with n a( j ~ 1, with the best fit n a( j = 0.96 (allow- 
ing for reionization would shift this even closer to 1). In 
large-field inflation models this translates into r < 0.3, so 
that the tensor contribution could be of the same order 
as the isocurvature contribution. Because of the degen- 
eracy between isocurvature and tensor perturbation such 
a situation is not altogether surprising. 

The COBE data together with as have an important 
role in constraining the spectral indices of the models. 
The Boomerang and Maxima-1 data are important in 
constraining the amplitude of the isocurvature contribu- 
tion, as they outline the pattern of acoustic peaks and 
thus are able to distinguish between the different patterns 
of peaks ]|l| in the isocurvature vs adiabatic models. 
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FIG. 4. The data points of COBE (o), Boomerang (•), and 
Maxima-1 (o), and the angular power spectra of five mod- 
els: (a) Our best-fit mixed (adiabatic+isocurvature) model 
(model 3, the solid line with the maximum at I ~ 200). (b) 
The best fit adiabatic model (model 1, dashed), (c) The 
best-fit isocurvature model (model 2, the solid line with the 
maximum at I ~ 300). (d) A mixed model with the largest 
(Ax 2 = 4) allowed a = 0.63 (model 5, dot-dashed), (e) A 
mixed model with the largest allowed «200 = 0.13 (model 6, 
dotted). 

V. CONCLUSIONS 

The main conclusion of the present work is that the 
Boomerang and Maxima-1 data definitely rule out all 
purely isocurvature models, including those with a large 
tilt |l7],[l8|]. This is because Boomerang and Maxima-1 
define the shape of the angular power spectrum in the 
region of the first acoustic peaks much more accurately 
than the previous data. 

However, a significant isocurvature contribution re- 
mains allowed. Indeed, as much as about half of the 
power at low multipoles could come from the isocurva- 
ture contribution, although the adiabatic contribution 
must dominate the first acoustic peak. This certainly 
leaves much room for various particle physics models. We 
should however point out that in many models, such as 
two-field inflation models or Affleck-Dine models, isocur- 
vature and adiabatic fluctuations are not completely in- 
dependent, and much more stringent constraints on a 
(and possibly nj so ) could be obtained in specific cases. 

Note also that we cannot obtain any clear conclusion 
about the spectral index n; so ; except for large n; so , the 
correlation between a and ni so is small. MAP and Planck 
will provide new accurate data on high multipoles and 
therefore could further constrain n- lS0l but at moderate 
tilts we do not expect MAP to improve much on the 
Boomerang and Maxima-1 limit on a. The situation 
with Planck is different, provided that high quality po- 
larization data is achieved as expected. Polarization will 
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FIG. 5. The angular power spectrum, Ci, and the adiabatic 
(dashed) and isocurvature (dot-dashed) contributions to it, 
Cf d and Cf°, for model 5 with (n ad ,n iso ) = (1.20,0.90) and 
a = 0.63. 
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FIG. 6. Same as Fig. [s| but for model 6 with 
(n ad ,n iso ) = (1.00,2.06) and a 2 oo = 0.13. 



resolve the degeneracy between isocurvature and tensor 
perturbations and push the limit on a lower by an or- 
der of magnitude (l9|. Nevertheless, as discussed in the 
present paper, already the balloon experiments on the 
temperature fluctuations of the microwave sky can yield 
interesting constraints on particle physics models. 
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